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Abstract 



Motivated by a recent tunneling experiment in a double quantum-well system, which reports 
an anomalously enhanced electronic scattering rate in a clean two-dimensional electron gas, 
we calculate the inelastic quasiparticle lifetime due to electron-electron interaction in a sin- 
gle loop dynamically screened Coulomb interaction within the random-phase-approximation. 
We obtain excellent quantitative agreement with the inelastic scattering rates in the tun- 
neling experiment without any adjustable parameter, finding that the reported large (> a 
factor of six) disagreement between theory and experiment arises from quantitative errors 
in the existing theoretical work and from the off-shell energy dependence of the electron 

self-energy. 
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A central quantity in the theory of interacting electron systems is the quasiparticle 
lifetime, which is the inverse of the scattering rate or the broadening of the quasiparticle 
state, and therefore, determines the width of the quasiparticle spectral function. The concept 
of inelastic lifetime is also important in electronic device operation because it controls the 
electron energy dissipation rate. It is, therefore, of great significance that a recent direct 
measurement of inelastic broadening in a two-dimensional electron gas by Murphy et al. 
reports a factor of six discrepancy between experimental results and the existing theory. In 
this letter, we develop a theory for inelastic Coulomb scattering lifetime in a degenerate 
two-dimensional electron system (2DES), finding essentially exact quantitative agreement 
with the tunneling results reported in Ref. 0. We also identify the reason for the factor of 
six disagreement reported in Ref. 

Over the past several decades two-dimensional electron systems have been extensively 
studied for both their fundamental and technological interest. The 2DES in high mobil- 
ity GaAs/AlGaAs heterostructures has become an especially suitable system for studying 
electron-electron interaction effects because of the reduced effect of impurity scattering aris- 
ing from the modulation-doping technique. Many properties of the 2DES are strongly influ- 
enced by the presence of electron-electron interactions. One important such property is the 
broadening of the electronic states by inelastic Coulomb scattering, which plays a major role 
in many physical processes, such as tunneling , ballistic hot electron effects , transport 
0, and localization The asymptotic properties of Coulomb scattering in a 2DES are 
well established from the existing theoretical work |l5HT0[: The electron inelastic lifetime 
Te in a pure 2DES becomes t~^(^) oc ^^In^ ioi Sp ^ i ^ ksT, and t~^(T) oc T^lnT for 
£f ^ ksT ^ ^, where ^ is the quasiparticle energy with respect to the Fermi energy ef, 

and T are the Boltzmann constant and temperature, respectively. Earlier experimental 
work on the inelastic lifetime of 2D electrons focused on the dephasing time |Tl|, while the 
recent experiment on tunneling in a double quantum-well structure directly measures the 
inelastic broadening. One advantage of the tunneling experiment over the dephasing exper- 
iment in this context is that the subtlety associated with quantum interference effects can 
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be avoided in a tunneling experiment, which directly obtains the inelastic broadening. In 
this sense the lifetime measured from the tunneling experiment is an excellent candidate for 
a direct comparison with theoretical calculations. It is the aim of this work to calculate the 
inelastic quasiparticle lifetime due to the Coulomb interaction in a clean 2DES and compare 
it with the results of the tunneling experiment 

The scattering rate obtained from the tunneling experiment with the contribution 
from the residual impurity scattering excluded, is essentially due to electron-electron inter- 
action. The effect of phonon scattering |]12|, including both acoustic and LO phonons, are 



safely negligible [||] in the experimental temperature range. Unexpectedly, (and as mentioned 
above), a very large qitative disagreement between the tunneling experiment and the existing 
theoretical calculations was reported. For example, the measured Coulomb scattering rate 
m close to the Fermi surface at low temperatures is found to be more than six times larger 
than that of the quoted calculation of Giuliani and Quinn (GQ) 0, which has been exten- 
sively used to interpret experimental results P, pn] , p!3| . This level of discrepancy is difficult to 



understand since the essential approximation used in GQ's calculation is the random phase 
approximation (RPA), and the corresponding three-dimensional RPA Coulomb scattering 
calculations are in excellent agreement with experiments |T^JT^. The biggest inaccuracy 
in the RPA comes from treating the short range correlations poorly. These short range 
correlations should not be very important in the low temperature scattering rate of elec- 
trons close to the Fermi surface where only long wavelength excitations are involved. This 
large discrepancy, if proved true, would cast serious doubt on the validity of the schemes of 
the existing theoretical work. This is particularly important in view of recent suggestions 
T6| that an interacting 2DES may not be a Fermi liquid and may have non-perturbative 



similar interaction effects akin to Luttinger liquids ||T^. Another significant puzzle is that 
the Coulomb scattering rate measured at very low temperatures as a function of energy 
in the quantum interference experiment ||Tl| seems to agree quite well with the GQ result. 
Since the calculation of GQ is the most widely used theoretical result in this subject, it is 
of considerable importance to investigate this discrepancy. For this purpose, we calculate 



the inelastic lifetime by obtaining imaginary part of the electron self-energy using the RPA 
dynamically screened exchange interaction, which is the same level of approximation as the 
work of GQ. Our calculation, with all the input parameters taken from the real samples, 
ie. with no adjustable parameters, shows very good quantitative agreement with the tun- 
neling experiment. We further find that the originally reported large discrepancy is due to 
some quantitative errors in the previous theoretical work and the negligence of the energy 
dependence of the Coulomb scattering rate. The Coulomb scattering broadening studied 
in this work is caused almost entirely by quasi-particle excitations. Plasmon excitation can 
contribute only at much higher electron energies 

We first present our calculation of the electron self-energy for a pure and ideal 2DES. 
The corrections from finite well-thickness, vertex correction, diffusive effects, and phonon 
scattering, all of which are included in our numerical work, will be briefly discussed at the 
end. The finite temperature electron self-energy in the RPA is 

ci P iuik 

where (3~^ = ksT, v is the area of the 2D system, and Kc = '^^(q') [^(q', '^t^)]^^ and Q° are 
respectively the screened Coulomb interaction and the electronic Green's function. After a 



standard procedure of analytical continuation |T5|,|T8[, the imaginary part of the self energy 
is obtained as 

1 1 
ImE(/c,^) = -^t;(g)Im— — — — [n^(^q+k) + nB(^q+k - w)], (1) 

^ q '?q+k - + ^U^j 

where = Yi'k^ j^m* — sp^ is the electron energy relative to the Fermi energy, v^p(b) is the 
fermion (boson) distribution function nj7(5)(x) = [e'^'^il]^^, v[(i) = 2Txe^ /tgq is the Coulomb 
potential. The RPA dielectric function is 

e{q,uj) = l-v{q)xl{q,uj), (2) 

where 

XM^^) = — —■ ] — 7 — (3) 



with 7 related to the mobihty broadening by 7 = e/{m*fi). The above particle-conserving 
polarizability [|1^] includes the essential effect of disorder scattering: The motion of electrons 
becomes diffusive rather than ballistic at large time and length scales. This expression allows 
a simple quantitative treatment of the diffusive effect arising from the finite value of mobility 
due to impurity and phonon scattering. For the experimental high mobility samples, the 
low temperature mobility is high, /i > lO^cm^/Vs, making the effect of disorder and phonon 
scattering practically negligible |[l|JT^. The use of Xc does not change the results within the 
numerical accuracy, however, it helps to improve the numerical integrations by suppressing 
the singularities associated with plasmon excitation. The noninteracting density-density 
response x° iii the above expression is 

Xo{q, ^ + 27 = - —-7 c , ■ ■ ^ 

From Eqs. (|l])-(^, ImE(A;, u;) can be computed. It is helpful to make clear the relation- 
ship between the lifetime obtained from this self-energy and the lifetimes calculated from 
the Fermi's Golden rule before we move on to discuss the numerical result. The lifetimes of 
electrons and holes from the Golden rule are 

2tt 1 

T"e"^(^) = y ^ nF<7'(^p)[l - nj.^/(^p_q)][l - n^^(^k+q)]|VkqP5(^k+q + -Cp-q " 6 " Q, 

2n 1 

rh\k) = —— ^ n^^/(^p)[l - nj.^,(^p_q)]n^^(^k+q)|VkqP5(^k+q - ^p-q + 6 " ^p), 
pqff' 

with Vkq = v{q)/e{q,C,k — ^k+q)- The last equation above is the equilibrium condition. 
Defining the broadening r{k,^k) = — 2ImS(/c, it is straightforward to show |18| 



It is therefore clear that the lifetime r~^{k,^k) = [— 2ImE(A;, is the relaxation time 
of the electron momentum occupation number Uko-. In general, it differs from either the 
electron lifetime or the hole lifetime. In particular, F"^ and differ by a factor of 2 at 



the Fermi surface. It is readily recognized that the hfe time obtained from the measured 
spectral function in a tunneling experiment is F"^, not t^.. This 0, we believe, is one source 
of error (by a factor of 2) in interpreting the experimental results of Ref. 

In Fig. ^ we show respectively the numerical results of F as functions of temperature 
T, energy ^, and electron density Ng. Several of the familiar features are easy to see from 
the figure: F(A;i;',0) oc T^lnT for small T, F(A;,^) oc ^^In^ for small ^ at T = 0, and 
T{kFi^) oc 1/Ns at small T. It is interesting to compare the numerical results in Fig. |l] to 
the analytical expressions obtained from the low T and small ^ asymptotic expansions of 
ImS in Eq. ([l|): 

r(A:, 6) = ^ = -^(^)^ In ^ for » k^T » a, (7) 

r(A;,a) = -^(^)2ln^ for ef » ik » ksT. (8) 

AtiTl Sf £f 

The above asymptotic expressions are consistent with our full numerical results (the 
inserts of Fig. [1|). It should be noticed that the prefactor in Eq. (0) is different from that 
of the work by GQ The cause of this difference, we believe, is that the corresponding 
expansion of GQ is incorrect by a missing factor of {tt/2)'^. This can partially account (by 
providing a factor of 2.5) for the fact that many studies have reported a Coulomb scattering 
rate significantly larger than the GQ prediction. Note that the asymptotic expression in 
energy, Eq. (H), is the same as that in GQ, explaining the puzzle why measurement of 
energy dependent scattering rate agrees with the GQ result ||TT| . 

Next we directly compare our calculation with the recent tunneling experiment in 
a double quantum well system. For the case of equal electron densities on each layer, the 
tunneling current as a function of the external bias potential V is sharply peaked at = 0. 
The resonance width Vhwhm, the bias potential at half maximum, is a measure of the 
quasiparticle lifetime. Under the condition (satisfied in Ref. [Q) that the Fermi energy is 
much larger than all the other energy scales involved, the resonance width hT^ff = Vhwhm 
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T,jf = ^T{kF, 0) + ^T{kF, Vhwhm). (9) 

It is important to notice that the finite bias potential Vhwhm introduces an OFF-SHELL 
energy dependence into the scattering rate. This kind of energy dependence is a direct 
consequence of simultaneous momentum and energy conservation [lll,pO| in the tunneling 
process. Taking the value of Vhwhm along with the values of all other sample parameters 
from Ref. we numerically calculate Fe// as a function of temperature T |^. In order to 
make a realistic comparison, we include the effects of finite well-thickness, vertex correction, 
and a finite value of mobility. The method to include these effects is discussed below. The 
calculated Fg// is shown as the solid line in Fig. ^ together with the experimental data 
and the theoretical result of GQ (dashed line) quoted from Ref. [0. One can see that the 
agreement of the present calculation with the experiment is excellent. The discrepancy 
between the GQ result and the experiment on the other hand is very large. This large 
discrepancy is due partly to the error (a factor of 7r^/2 ~ 5) in GQ's work which was 
discussed following Eqs. (|^) and (|^, and partly to the negligence of the off-shell energy 
dependence of the scattering rate, which contributes a 30 ~ 40% quantitative effect. The 
excellent agreement between the present calculation and the experiment suggests that the 
commonly adapted Fermi liquid RPA-like many-body treatments for the Coulomb scattering 
rate are well valid in GaAs-based 2DES. This also shows that a clean interacting 2DES is, 
in fact, a Fermi liquid fl^ similar to a three dimensional system. 



Finally, we briefly discuss how the corrections from flnite well-thickness, vertex correc- 
tion, phonon and impurity scattering are incorporated into our calculation. The flnite well- 
thickness, which tends to weaken the interaction at short distances, can be represented by re- 
placing v{q) by v{q)F{q). The form factor may be chosen as = (2/g6)[l + l/g6(e-«^- 1)] 



with h as the well thickness ^2[. The influence of vertex correction, which tends to decrease 
screening through a local fleld correction, may be estimated by the replacement of x°c{li ^) by 
Xc(q', — G{q)] with the Hubbard local fleld approximation [|23l G{q) = 0.5q/ (g^ + kp)^/"^. 
The effect of LO phonon-mediated electron-electron interaction can be included [l^ by 
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adding the factor (1 — eoo/eo)/(eoo/eo — ^^^/^io) the dielectric function e{q,u!) with the 
parameters for GaAs materials as cq = 12.9, eoo = 10.9, and uj^o = 36.8 meV. For low 
energy excitations, the LO phonons act as a small source of static screening. The effect 
of a finite mobility due to impurity and acoustic phonon scattering, which is small under 
the present conditions, can be taken into account by putting into Xc appropriate value 
of 7 = e/m*fi. For computational reasons, all the data presented here are calculated with 
fi = lO^cmVV^s. 

In summary, we have calculated the inelastic scattering rate due to electron-electron 
interaction for a two-dimensional electron gas. Our work is motivated by the large disagree- 
ment between the recent tunneling experiment [|l| and the existing theoretical calculations. 
Using experimental sample parameters, we have obtained excellent quantitative agreement 
with the tunneling experiment. Our work suggests that a clean interacting 2DES is a Fermi 
liquid and that RPA-based perturbative many-body calculations are of quantitative validity. 

The authors thank J. P. Eisenstein for helpful discussions. This work is supported by the 
U.S.-ARO and the U.S.-ONR. 
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FIGURES 

FIG. 1. (a) The Coulomb scattering rate r{k,u) as a function of temperature T/Tp, where 
Tp = Ep/kB- The soUd hne is for T{kp,{)). The dashed hne is for T{k^^k) with = 0.46^. The 
insert shows A = hV{kp,{))£p /[{kBTf\^{Tp /T)] as a function of T/Tp. The density of the 2DES 
is Ng = 1.6 X lO^-'^cm"^. (b) The Coulomb scattering rate T{k,u!) as a function of electron energy 
^k/^F- The solid line is for r(fc,^fc) at T = 0. The dashed line is for T(k,$^i.) with T = O.lTp. The 
insert shows a = hr{k,^k)£F/[^k^'n-i^F/Ck)] as a function of Ck/^F at T = 0. The density of the 
2DES is Ns = 1.15 x lO^^cm-^. (c) T{kp,0)/T^ as a function of the density A^^ at T = 3K. 

FIG. 2. The Coulomb scattering determined tunneling resonance width Fe// as a function 
of temperature T. The solid hne is the present calculation with density Ng = 1.6 x lO^-^cm"^, 
well-thickness b = 200^", and mobility jj, = lO^cm"^ /Vs. The diamonds and the dashed line are 
respectively the experimental data and the theoretical result of GQ quoted from Ref. 1. 
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Fig.l(a) 
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Fig. 1(c) 
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Fig.2 
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